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Abstract. We analyze galactic black hole mergers and their emitted 
gravitational waves. Such mergers have typically unequal masses with mass ratio 
of the order 1/10. The emitted gravitational waves carry the imprint of spins 
and mass quadrupoles of the binary components. Among these contributions, 
we consider here the quasi-precessional evolution of the spins. A method of 
taking into account these third post-Newtonian (3PN) effects by renormalizing 
(redefining) the 1.5 PN and 2PN accurate spin contributions to the accumulated 
orbital phase is developed. 
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1. Introduction 

Mankind has admired the sky with its Milky Way, stars, planets, moons and the 
Sun, the meteors and the Northern Light. In the age of satellites we map the sky 
through all electromagnetic frequencies and also through high energy cosmic particles. 
Gravitational waves are the last great frontier to be reached and passed for the search 
for our physical understanding of the universe via the messengers reaching us with the 
speed of light or very close to it: 

Gravitational waves, these fleeting distortions of space, the signature of 
accelerating masses in the universe, remain to be discovered. Therefore it is of 
paramount importance to understand, which objects are about to generate strong 
gravitational waves. 

When the black holes residing in the centers of almost all galaxies merge, following 
the merger of their host galaxies, relatively strong gravitation waves are emitted and 
various aspects of this process were already investigated [l]-[6]. A recent review of the 
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generic aspects of these galaxies nuclei as sources for ultra high energy cosmic rays is 

Starting from the broken powerlaw mass distribution of galactic central black 
holes ISj-IHj, with simple assumptions on the cross section and relative velocity of 
the two galaxies, we have shown [6] that the most likely mass ratio v = m2/mi < 1 
of the merging black holes is between 1/3 and 1/30, thus a typical value will be at 
about 1/10. Within this range of mass ratios the key stages of the evolution leading to 
and following the merger can all be treated with post-Newtonian (PN) approximations. 
A numerical analysis of such mass ratios was recently presented in [12 ] . 

In the next section we describe the sequence of events which will accompany the 
merger of two supermassive galactic black holes. Then from Section 3 we concentrate 
on the gravitational radiation dominated dissipative regime. An important quantity 
in gravitational wave detection is the number of cycles spent in the frequency range 
of the detector. This quantity is proportional to the accumulated orbital phase (the 
total angle covered by the quasi-circular orbital evolution). The gravitational waves 
emerging from this process carry the imprint of the physical characteristics (finite size) 
of the black holes, like their spins and quadrupole moments. Spin precessions allow 
for increased accuracy and degeneracy breaking among the parameters characterizing 
the source [13]. 

We remark that according to [13] , the mass quadrupole - mass monopole coupling 
for maximally rotating black holes (with quadrupole moment originating entirely in 
their rotation) gives rise to contributions in the accumulated orbital phase larger by a 
half to one order of magnitude than the spin-spin contribution, at least for the typical 
unequal mass ratios. Therefore for a successful detection this qualifies as the second 
most important finite mass contribution to be considered, after the spin-orbit effect. 

In Section 3 we discuss the effect of the quasi-precessional evolution of the 
spins on the accumulated phase of the gravitational waves and we present estimates 
of the PN order in which various finite size contributions to the change in the 
relative configurations of the spins and orbital angular momentum contribute. The 
cumbersome equations leading to these estimates are presented in the Appendix. 

Sections 4 and 5 explore a convenient analytical way to encompass the effect of 
spin and quadrupolar quasi-precessions even at a lower level accuracy. The leading 
order contribution to the spin-spin (Section 4) and spin-orbit (Section 5) contributions 
to the accumulated orbital phase are discussed both for the equal mass and unequal 
mass cases in distinct subsections. 

We summarize our findings in the Concluding Remarks. 

2. The merger of two supermassive galactic black holes: an odyssey 

First the dance: Two galaxies with central black holes approach each other to within a 
distance where dynamical friction keeps them bound, spiraling into each other. In the 
case when both galaxies with their central black holes were radio galaxies, their jets 
get distorted and form the Z-shape [15]: the radio emitting blobs and tails produce 
the appearance of dancing veils. 

Second the meeting of the eyes: The central regions in each galaxy begin to act 
as one unit, in a sea of stars and dark matter of the other galaxy. During this phase, 
the central region can be stirred up, and produce a nuclear starburst. 

Third the lock: The black holes begin to lose orbital angular momentum due 
to the interaction with the nearby stars [16], and next by gravitational radiation. 
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The spin axes tumble and precess. This phase can be identified with the apparent 
superdisk, as the rapidly precessing jet produces the hydrodynamic equivalent of a 
powerful wind, by entraining the ambient hot gas, pushing the two radio lobes apart 
and so giving rise the a broad separation 17J. Due to the combined effect of precessions 
and orbital angular momentum dissipation by gravitational waves the spin (and hence 
jet) direction of the dominant black hole is reoriented approximately in the direction 
of the original orbital angular momentum, leading to a spin-flip [6 . 

Fourth the plunge: The two black holes actually merge, with not much angular 
orbital momentum left to be radiated away, thus the direction of the spin remained 
basically unchanged (in the case of the typical mass ratio range we discuss; for other 
setups the final spin may vary [18]-[23j). The final stage in this merger leads to a 
rapid increase in the frequency of the gravitational waves, called "chirping", but this 
chirping will depend on the angles involved. 

Fifth the rejuvenation: Now the newly oriented more massive merged black hole 
starts its accretion disk and jet anew, boring a new hole for the jets through its 
environment (the GigaHertz peaked sources). In this stage the newly active jet is 
boring through the interstellar clouds with a gigantic system of Shockwaves, that 
accelerate protons and other charged particles in a relativistic tennis-game. These very 
energetic particles then interact in those same clouds, that slow down the progression of 
the relativistic jet, and can so give rise to an abundance of neutrinos, TeV photons from 
pion decay, and yet other particles and photons. These particles and their interactions 
reach energies far beyond any Earth-bound accelerator. 

Sixth the phoenix: The newly oriented jets begin to show up over some kpc, and 
this corresponds to the X-shaped radio galaxies, while the old jets are fading but still 
visible [23]. This also explains many of the compact steep spectrum sources, with 
disjoint directions for the inner and outer jets. 

Seventh the cocoon: The old jets have faded, and are at most visible in the low 
radio frequency bubbly structures, such as seen for the Virgo cluster region around 
M87. The feeding is slowing down, while a powerful jet is still there, fed from the spin 
of the black hole. 

Eighth the wait: The feeding of the black hole is down to catching some gas out 
of a common red giant star wind as presumably is happening in our Galactic center. 
This stage seems to exist for all black holes, even at very low levels of activity. This 
black hole and its galaxy are ready for the next merger, with the next black hole 
suitor, which may or may not come. 

We predict that the super-disk radio galaxies should have large outer distortions 
in their radio images, and should be visible in low keV X-rays. A detection of the 
precessing jet, and its sudden weakening, would then immediately precede the actual 
merger of the two black holes, and so may be a predictor of the gravitational wave 
signal. 

3. The contribution of the spin precessions to the accumulated orbital 
phase 

In the rest of the paper we concentrate on the gravitational radiation dominated 
dissipative regime, when the supermassive black hole binary radiates away 
gravitational waves. 

If the binary system evolves unperturbed, the orbit circularizes faster than it 
shrinks in radius. Therefore we consider circular orbits for which the wave frequency 
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is twice the orbital frequency. 

Besides the Keplerian contribution to the accumulated orbital phase there 
are general relativistic contributions, which scale with the post-Newtonian (PN) 
parameter e = m/r ~ f ^ (with m the total mass and r the radius of the reduced mass 
orbit about the total mass; we take G — I — c). It is generally agreed from comparison 
with numerical evolution, that all contributions up to 3.5 PN orders have to be taken 
into account. To this order, beside the general relativistic contributions starting 
at IPN, there are various contributions originating in the finite size of the binary 
components. These are related either to their rotation (spin) or to their irregularities 
in shape (like a mass quadrupole moment). 

Due to the continued accretion of surrounding matter into the supermassive 
galactic black holes they spin up considerably. By taking into account only the 
angular momentum transfer from accreting matter, the dimensionless spin parameters 
Xi = Si /ml (with i — 1,2) grow to the maximally allowed value 1 even for initially 
non-rotating black holes [55]. By taking into account the torque produced by the 
energy input of the in-falling (horizon-crossing) photons emitted from the steady-state 
thin accretion disk, which counteracts the torque due to mass accretion, the limiting 
value of the dimensionless spin parameter is slightly reduced to 0.9982 [26]. Various 
refinements of this process, with the inclusion of open or closed magnetic field lines |27j - 
[55] . also jets in the magnetosphere of the hole [35]- [37] have not changed essentially 
this prediction. Therefore we assume for the present considerations maximal rotation 
Xi ss 1. Due to this rotation the black hole is centrifugally flattened (it becomes an 
oblate spheroid), a deformation which can be characterized by a mass quadrupole 
scalar Qi [38], or its dimensionless counterpart pi — Qi/mirn^ « — (nii/m) . 

The accumulated orbital phase can be formally given as the PN expansion [T4] . 
[39], [40] 

(j) = (l)c + 4>N + <PlPN + <Pl.5PN + 4>2PN + 4>2.bPN + 03PAf + <p3.5PN •(!) 

Spin-orbit (SO) contributions appear at 1.5 PN, spin-spin (SS; composed of proper 
S1S2 and self-contributions Sf and 83) and mass quadrupole - mass monopole coupling 
(QM) contributions at 2PN. They come with various PN corrections (PNSO at 2.5 
PN; PNSS, PNQM and SO^ at 3PN; finally 2PNS0 at 3.5 PN; these were not 
computed yet with the exception of PNSO [1I]-[12|). At 3PN there are additional 
spin and quadrupole contributions, originating in the quasi-precessional evolution of 
the spin and orbital angular momentum vectors. The number of gravitational wave 
cycles J\f can be computed from the accumulated orbital phase as A/" = (0c — 4>) 
(the gravitational wave frequency being twice the orbital frequency for quasi-circular 
motion). 

The spin precession equations with SO, SS and QM contributions were given by 
Barker and O'Connell [43]. The relative geometry of the two spins Si and orbital 
angular momentum L can be best described by the angles 7 = arccos(Si • S2) and 
Ki = arccos(Si • Ln) (an overhat denotes the direction of the respective vector). These 
angles are related by the the spherical cosine identity 

cos 7 = cos Ki cos K2 + cos Alp sin ki sin K2 , (2) 

with Alp the relative azimuthal angle of the spins. Due to quasi-precessions, Hi and 7 
evolve. Supplementing the SO and SS contributions already given [44] with the QM 
contributions (with the quadrupole moment arising from pure rotation), we give these 
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expressions to second order accuracy: 

(cOSKi) = (cOSKj)sO + {cos Ki)'ss + (cOSKi)^Af 

(cos 7)" = (cos7)so + (cos 7)55 + (cos7)qa/ > 



(3) 
(4) 



with the detailed contributions enhsted in the Appendix. 

The order of magnitude estimates of the terms in Eqs. (jA.l[) - (|A.9p (cf. a footnote 
in Ref. [45], according to which 0{dx) = 0{x)re^^^^), and assuming maximal rotation 
lead to 



Thus the angle 7 varies at IPN, while the angles Kj only at 1.5 PN during one 
orbital period. However 7 appears only in the 2PN order S1S2 contribution to the 
accumulated orbital phase as opposed to the angles Ki which are present in all finite- 
size contributions, in particular in the 1.5 PN order SO contribution. 

Therefore 3PN order quasi-precessional contributions to the accumulated orbital 
phase in the comparable mass case v 1 arise from Ssoi^i, Sss^i, SqmHi and Ssol] 
while in the unequal mass case, found typical for supermassive galactic black hole 
binaries only from 630^2, ^33^2 and Ssol- 

4. Leading order evolution of the spin-spin coefficient 

From among the quasi-precessional evolutions, we have recently integrated [14j the 
SO evolution (averaged over one orbit) of the angle Aip, finding 



Here /i = mim2/m the reduced mass, a the radius and n = 21: /TorUt = jT^ave- The 
time-dependent expression for 7 is then given by Eq. ([2]), with A"?/; from Eq. ([6]). The 
variation of both A?/) and 7 thus have a periodicity with period 
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Figure 1. (Color online) (a) The factor F (y), shown for the whole allowed range 
of mass ratios, diverges in the equal mass case, (b) The factor F (y) in the typical 
mass ratio range (1/30, 1/3) found for merging supermassive galactic black hole 
binaries is of order unity. 

Equal masses 

For equal mass binaries this time-scale becomes infinite liuii^^i F (i/) — (x (see Figla) 
expressing the fact that the time-dependent part of 7 goes to zero. Thus 7 is a constant 
and there is no 3PN quasi-precessional contribution to the accumulated orbital phase 
for equal mass binaries. 



4-2. Unequal masses 

For the typical mass ratio the factor F (i^) is of the order of unity (see Fig lb), thus 
the variation time-scale of 7 is e^^ G (10,1000) times larger than the wave period. 
The contribution of the quasi-precessions can be taken into account then by keeping 
only the first term in Eq. Q when computing the 2PN order S1S2 contribution to 
the accumulated orbital phase, e. g. by introducing a renormalized spin coefficient. 
This coefficient arises by modifying the 2PN contribution with the average of some 
of the 3PN contributions. By using this renormalized coefficient in the context of the 
2PN accurate dynamics, one can bring closer the 2PN analytical prediction to the 
numerical results 14jl3 

Thus, in the unequal mass case we propose to replace 

'S'iS'2 
4877TO'' 



o'SiS2 = 73 j(— 247cos7 4- 721 cosKi cos K2) , (9) 



with 



'^SiS2 = o^^4 "^"^ ^1 ^2 , (10) 



79S'iS'2 
Srjm 

in the respective 2PN contribution to the accumulated orbital phase 

1 / 9275495 284875 1855 , 15 \ 1/8 

<t>2PN = 1 V H fl 0- T^'^ , 

T] V 14450688 258048 ' 2048 ' 64 y 

cr = CrSiS2 + (^SS-self + <yQM ■ (H) 

(Here 77 = /i/m is the symmetric mass ratio and r a dimensionless time parameter.) 



I However we stress that there are other finite size 3PN order contributions to the phase, like PNSS, 
PNQM and SO^, which are still not computed. 
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5. Leading order evolution of the spin-orbit coefficient 



-/3-7r 



-1/4 



The 1.5 PN contribution to the accumulated orbital phase is 
with the SO contribution given by 



(12) 



(13) 



5.1. Equal masses 

The quasi-precessional equations with the SO, SS and QM contributions included were 
recently integrated by Racine for the equal mass case According to his analysis, 
there is a previously unknown constant of the motion 
2 

(^i cosKi + 5*2 C0SK2) ■ (14) 



A 



L 



N 



As for equal masses P = {AT L n / 24) A =:const., the quasi-precessional evolution at 3PN 
does not affect the coefficient (3 up to the 3PN accuracy. 



5.2. Unequal masses 

In the unequal mass case the leading order evolutions of the angles Ki are 8301^2 and 
5 33^2- Averaged over one orbit they give [44] 



K2 



2a3 



■ sin Ki sin Aij) 



(15) 



For unequal masses v k, 10^^ the estimates ([S]) allow to consider ki as a constant (in 
a first approximation) . The leading-order time-dependence of Aip is given by Eq. ([6]) . 
Integration gives 



'«2 



('^2)0 



.1/2 



sm Ki cos 



(AV')o 



2a 



(16) 



(We have employed e ~ m/a .) The time-varying part has the same period T3PN3 
given by Eq. ([7]), which is of the scale of e^^ the orbital period. In agreement with 
this and the estimate ^ 



0('*2-(K2)o)_eV2 



(17) 



O {630,331^2) £3/2 

As the time varying part averages out, to leading order accuracy we can simply identify 
the angle K2 with its constant initial value (1^2) q. 

In conclusion, in the unequal mass case and in the leading order approximation 
Ki =const. the time variation due to quasi-precessions renders K2 to another constant 
at 3PN accuracy. There is no need to renormalize (3, when taking into account the 
leading order quasi-precessional contributions. 
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6. Concluding Remarks 

The process of the merger of two supermassive galactic black holes gives rise to exciting 
signals in all the messengers accessible now and in the near future, and will allow 
to probe some of the most profound secrets of Nature, occurring at energies and 
circumstances far beyond anything on Earth. We have described the sequence of events 
during the process of merging of such supermassive galactic black holes with reference 
to related processes involving accretion, jets, energetic particles and electromagnetic 
signatures. We argued that the class of radio galaxies with a super-disk are particularly 
promising candidates for future powerful low frequency gravitational wave emission. 

The accumulated orbital phase is a quantity to be known precisely for successful 
detection of gravitational waves. Besides Keplerian, first and second PN order general 
relativistic contributions, the finite size of the black holes also contributes, first at 
1.5 PN orders via the spin-orbit interaction, then at 2PN by the spin-spin and 
quadrupole-monopole coupling. These contributions are further modulated by quasi- 
precessional evolutions of the spins. However this is but a subset of all possible finite 
size contributions at 3PN, as stressed earlier. The way we propose to take into account 
the precessional 3PN contributions is to add their average to the lower order terms, a 
procedure we call renormalization of the lower order coefficients. 

We have proven here that there are no such quasi-precessional contributions in 
the equal mass case. Furthermore we have proved for the unequal mass case that 
renormalizing the 1.5 PN spin-orbit coefficient will not change its value. At 2PN 
there are quadrupolar and spin-spin contributions to consider. As the former can be 
regarded as constant up to 3PN, the only modification due to quasi-precessions comes 
from the spin-spin interaction. The renormalized spin-spin coefficient at 2PN thus is 
adjusted by the average of the 3PN quasi-precessional contribution, taken over the 
period of change induced by these quasi-precessions (a time-scale one PN order lower 
than the orbital period). We have checked that the renormalized spin-spin coefficient 
gives closer results to the numerical estimates for compact binaries with the mass ratio 
of order 1/10 [I^. 

Our results related to the renormalization of the spin coefficients are significant 
for supermassive black hole binaries, which are typically LISA sources. However as we 
have employed nothing but the mass ratio of the binary in our considerations, they 
also apply for binaries composed of an intermediate mass and a stellar mass black 
hole, which qualify as LIGO sources. 
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Appendix A. Quasi-precessional evolutions of the relative spin angles 



The SO and SS contributions of the instantaneous variations of the relative spin angles 
[44] supplemented with the respective QM contributions (assuming the quadrupole 
moment arises from pure rotation), to second order accuracy are: 



(cOSKi)so 

(cos K,i)'ss 
(coski)qa/ 



(cOSK2)so 
(cOSK2)ss 

(cOSK2)qm 

(cos7)so = 
(cos7)ss = 



(2 + Ln • (Si X S2) , 
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(f •S2)Ln- f xSi 



^(f -Si)? .(Si XS2) 



P2 



|i(f •Si)LN-(f xSi) 



L 



N 



(f • S2)f • (Si X S2) 



0(2 + i.)LN-(Si XS2) 
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(f •Si)LN-(f XS2) 



(f •S2)f -(Si XS2) 



^(f •S2)LN-(f XS2) 

J2 



+ |^(f -Si)? • (Si X S2) 



2r3 
3 



(z.-i/-i)Ln-(Si XS2) , 
52(f ■S2)-5i(f -Si)] f.(Si X S2) 



(cos 7); 



IQM 



^(r-S2)-— (r-Si) 



f-(Si X S2) 



(A.l) 
(A.2) 

(A.3) 
(A.4) 

(A.5) 

(A.6) 

(A.7) 
(A.8) 
(A.9) 
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